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CHORDWISE AND COMPRESSIBILITY CORRECTIONS TO SLENDER-WING THEORY ! 


By Harvarp Lowax and Loma SrnupER 


SUMMARY 


Corrections to the solutions given by slender-wing theory for 
the lift distribution on triangular and rectangular wings of low 
aspect ratio are obtained in two steps. First, slender wing 
theory 18 used to find the load distribution over a wing of giten 
shape. Second, the spanwise variation of the loading so ob- 
tained is left unchanged but the chordwise variation is modified 
by satisfying an appropriate integral equation. Results are 
shown for flat-plate wings and, in the case of ihe subsonic, 
iriangular wing, а comparison is made with other theoretical 
solutions and experimental results. 


INTRODUCTION 


The calculation of loading on three-dimensional lifting 


surfaces is a fundamental problem in aerodynamic research. 
The complexity of the problem has led to the development 
of certain simplified theories by means of which the loading 
on special types of plan forms can be estimated quickly. 
The amount of error which these estimates contain is of 
considerable interest, as are methods which will tend to 
correct, such errors without undue labor. 

Slender-wing theory applies to one, such simplified body 
of analysis. "There are two basic assumptions of this theory. 
One, the angle of attaek is small enough so that the vortex 
sheet lies in the plane of the wing and the boundary condi- 
tions for the wing can be projected onto & horizontal plane 
parallel to the direction of the free stream; and the other, 
that either the chord wise gradient of velocity is small enough 
or the free-stream Mach number is close enough to unity 
that the linearized partial differential equation which governs 
the fluid flow becomes Laplace’s equation in a plane trans- 
verse to the free-stream direction. References 1 through 6 
are examples of papers developing slender-wing theory. 

The magnitude of the error of such & theory, in the case 
of subsonic flow, is indicated by observing solutions for 
triangular wings. Slender-wing theory gives a finite value 
for the loading along the trailing edge. Proper inclusion of 
the chordwise and compressibility effects results in solutions 
that satisfy the Kutte condition and make the loading fall 
to zero at the trailing edge. It is the purpose of this report 
to study such modifications. 

The corrections due to the chordwise and compressibility 
effects are obtained in the following manner: First, an in- 


t Supersedes NACA TN 2295, “Chordwise and Compressiblity Corrections to Slender- 
Wing Theory," by Harvard Lomar and Loma &luder, 1961. 
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tegral equation is set up relating the shape of the wing 
surface to the lift distribution; second, this integral equation 
is solved for the given shape under the assumption that the 
chordwise velocity gradients are small or that the free-stream 
Mach number is unity; and finally, the integral equation is 
reinspected for the same wing shape, this time with the 
spanwise lift distribution fixed at the variation just obtained 
&nd with the chordwise variation as the unknown and the 
Mach number terms included. | 

Results are presented and discussed both for trianguler | 
and rectangular, fl&i-pl&te plan forms in both subsonie and 
supersonio flow. | 


LIST OF IMPORTANT SYMBOLS 


. fb 
4 &spect ratio (5) 
b span of wing measured normal to plane of 
symmetry 
Co root chord of wing 
Cz lift coefficient 28 
C pitching-moment coefficient 
pitching moment about leading edge or apex 
@абев 
E(t,k) incomplete elliptic integral of the second kind 


with argument ¢ and modulus Ё 


| Be p= [+ E E 


complete Ma Seo of second kind, Е (1,4) 
incomplete elliptic integral of the first kind 
with argument # modulus Ё 


Gk 

| '«5- [Sam | 

K complete elliptic integral of first kind, F(1,k) 
m 


for triangular wing, slope of leading edge relative 


to plane of symmetry 
Af, free-stream Mach number 





E 
F(tk) 


Vo 
(seed of sound in free stream 


p static pressure 

Àp Piı— Pu 

q free-stream dynamic pressure (5 po V. 2) 
8 Semispan of rectangular wing 

S area of wing 
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u perturbation velocity component in the direction 
of the x axis 
Au Uy — i 
Vo free-stream velocity 
w perturbation velocity component in the direction 
of the z axis 
Wo —Voa 
г,у,2 Cartesian coordinates of an arbitrary point 
Zu Cartesian coordinates of source or doublet 
position 
£ 
To Co 
a &ngle of attack 
8 4 —AMg[ 
05 pm 
po density in free stream 
T, Фр doublet weighting factors 
e perturbation velocity potential 
Ae Pu— Фі 
SUBSCRIPTS 
i conditions on lower surface of wing (at 2==0—) 
ч conditions on upper surface of wing (at z=0 +) 


THE INTEGRAL EQUATIONS 
SUBSONIC 


Triangular plan form.—<A general solution of Laplace's 
equation which is suited to problems in linearized subsonic 
wing theory (given, e. g., in reference 7) is that which relates 
& velocity potentiel or perturbation velocity to the value of 
its jump across а given surface. For the lifting triangular 
wing shown in figure 1 this can be written 


Au dY, 
=) da Bog zi 8Xy— VY BEER Q 


where 8—4/11— Mil, v is the perturbation velocity parallel 


to the x axis, and Au is the jump in е over the wing plan form. 





FIGURE 1,—Triangular-wing coordinate system. 
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In linearized theory this jump сап be related to the loading 
coefficient Apg by the equation 


20) . K 


Further, the velocity potential e can be found by the relation 


e- | мат 


Operating on equation (1) in this manner and interchanging 
the order of integration gives 


тї dyi 
eras fo an 


42 o. T 3 
(т-ту e+ By — yo 3) 


which represents, physically, a distribution of elementary 
horseshoe vortices. 

The effect of compressibility in a lincarized study of lifting- 
surface theory can only enter through the use of 8. Setting 


T—H 
TUS | " фа 
~~ = — tif 840 — y)? -- B72? x 


it is seen that с is the only term in equation (3) which con- 
tains B. This term has an Interesting interpretation in the 
light of the study which has been made at sonic speeds. At 
Mal (i. e., 8—0), с has either the value 2 or 0, depending 
on whether x, is less or greater than z. Hence, for Af,—1, 
equation (3) becomes 


MEL s [7 an P\ ay, 


gu er 5) 


Now reversing the order of integration and using the define- 
ton implied by equation (2), namely, 


2р 2 дАе 
Fo Әд 
gives finally 
- |" Ag di | (6) 
97 2x -my =NH г* 
Equation (6) has been studied in reference 5 as the funda- 
mental equation for slender wings or wings flying at near 
sonic speeds. It is an equation which gives the solution for 
the velocity potential in a three-dimensional flow in terms of 
two-dimensional doublets, the two dimensions being at right. 
angles to the free-stream direction. A solution of such a 
nature is immediately implied by the physical character of 
both sonic wing theory, in which the Mach cone has degener- 
ated to a Mach plane, and slender-wing theory, in which the 
wing is so slender that the chordwise gradient of velocities 
can be neglected compared {о Ше vertical and lateral 
gradients. 

By comparison of equation (5) with equation (3), it is 
seen that the term e can be interpreted as & factor which 
corrects the slender-wing-theory results as given by equation 
(5) for the effects of chordwise gradients in velocity and 
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Fiat RE 2.—Typical тагіа ор ofe over a triangular wing. 


compressibility. By consideration of the effect at one point 
of the distribution of doublets over the wing, this correction 
ean be visualized as & reweighting of the two-dimensional 
doublets according to their position relative to the point. 
Figure 2 indicates the variation of с across the span at various 
chord stations for 8—0.6. Observe that the doublets ahead 
of the point at which the potential is to be determined are 
still weighted far more heavily than those behind the point. 
The effect of considering В different from zero, however, is 
to reduce the extreme difference in weight occasioned at 
8=0 so that the doublets behind a given point do have some 
effect on the induced velocities there, and the doublets ahead 
of a point induce a somewhat smaller disturbance than 
before. Since the strength of these weighted two-dimensional 
doublets is given by the magnitude of the three-dimensional 
loading, their strength is zero everywhere off the wing plan 
form including the area behind the wing occupied by the 
vortex wake. 

Two different methods for the further reduction of equation 
(3) will be considered. The first method involves finding the 
vertical induced velocity for points along the z axis, while 
the second involves finding the average vertical induced 
velocity along the span at a given chord station. The first 
method must be discarded for triangular wings because of 
difficulties around the apex; the second, however, proves to 
be satisfactory. The simplification obtained by considering 
the vertical induced velocity for points along the x axis will 
be considered later in connection with the rectangular wing. 

Since it is easier to consider first the averaging process, 
the operator 

lim. 1 9 ("* 
24 IMT DZ lome 


is applied to the weighted doublets, ez/[(y—y:) +27, of equa- 
tion (3} witk the result that БО 


I 
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Ар Ар 
rece, pan, (2) an (7) 
sel in V mi—y 8x Jo -mz 2mx(z—2) 


Еу р | уши meray A 
mr— yı mz-l-yi 


where 10 is the average value of the vertical induced velocity _ 

along a given span. | NT 
The solution for Ap/g obtained from slender-wing theory 

сап be written ? | 


Vo 


w= 








са 
di; 


Ap 4wo mi 2; 
та ^ (Е 
where in that theory №. (2/со)=1. И the value of Ap/¢ 
given by equation (8) is placed in equation (7), the resulting 
integral equation can be written in a simplified form if it is 
noted that l 


(Jona т — 01 mz4-yi 


(8) 


ү Ү(ж—х)*--8°'4т 
отет _.. 
ma-r) ү — m'(z!— у )-Е2ттх — т? 


where for the first term in the brackets the transformation 


(9) 


7;—mzr—j, was used and for the second the transformation — 


=т=. Hence, equation (7) finally reduces to the 
following 


E | [ AR) ize (* zihÁ (2) ч (10) 


2тх 4 T?— r? (z—21) 


The solution'of equation (10) will be discussed in a later 

section devoted to triangular wings. E 

Rectangular plan form.—If the plan form of the wing is 

rectangular as shown in figure 3, then equation (3) is modified 
slightly to the form 

. A 
JA de 
9 8r Jo yy Tz 


It is possible in this case to study the vertical induced veloc- - 
ity for points along the хт axis; that is, to find 24/92 by 
equation (11) and then set both y and z equal to zero. In 
order to do this & special notation is employed. Thus, if 
the indefinite integral of f(y)/y* can be written (where f(y) is 
bounded at y—0) 


(11) 


[9 ay- sec 
then, by definition, 
+ ae dy=J(e)—J(—s) (12) 
By means of this definition it can be shown (see reference 8) 


that 
* This solution follows from an analysis of equation (6). Bce reference 5. 
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peres 


K 


FiavRES3.—HRoctangular-wing coordinate system .- 


wma |, 7 aan Pel кш zeros 


If? 


АЙ шы 


2 (14) 


Ao СА f. 
Voto һ» (=) y 


integration by parts gives 


82—12 
Ü TB yi 


=) а 


—c 


E 


and equation (13) can be written 


waza |“ (= aum з) Z) das, 


This integral equation has been derived previously by K. 
Wieghardt (reference 9) with regard to the rectangular-wing 
problem. The solution of equation (16) will be discussed 
in & later section devoted to rectangular wings. 


| н ә 


(16) 


SUPERSONIC 


Triangular plan form.—In passing from subsonic to 
supersonic theory, we pass from the elliptic to the hyperbolic 
partial differential equation and in particular fram Laplace’s 


1 The solution for the rectangular wing given by slender-wing theory is that the load be 
zero across every spanwise strip aft of the leading edge. To find the chordwise correction to 
such a theory, therefore, a spanwise distribution must be assumed, Since, however, slender- 
wing theory also requires an elliptical span loading for the boundary conditions of a rectan- 
gular wing to be satisfied, a reasonable choice is that given by equation (14). 
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equation to the wave equation. The solution which relates 
the perturbation velocity u at any point in the field to the 
loading on the wing can again be written in terms of an 
elementary. horseshoe vortex distribution over the wing plan 
form. As in reference 8, this becomes 


2 = |, (r— —t)Audzidy: rr 
— 2v Ot JJ. куу) zx —2) — By — "2 pts? 





Ap _2Au 
gq Vo 





and since u—2* an d — 
T 


Ap dzidyi 


Voz [Г gro | f—1 
Tagan yg 17 
um (y— -yy + 2? Кх — — 1)? — a B'ty — у). 2 ( 
where т is the area on the wing bounded by Ње edges and 
the trace of the Mach forecone from the point x, y, 2. Again 
the effect of compressibility appears only in the term within 
the brackets. Hence, setting 


— &— 7 ' 
” AG—ny—8Hy—yy—gizi 


с, contains all of the Mach number effects at supersonic 
speeds. At M,—1, e,—1, and since, by the definition of 
т, 21 €z, it follows that at sonic speeds equation (17) also 
reduces to equation (5). Therefore the doublet distribu- 
tions represented by equations (17) and (3) аге consistent 
at the speed of sound. 

In order that an exact parallel can be provided with the 
subsonic solution to the triangular wing, the average vertical 
induced velocity for points along a given span is again 
considered. It сар be shown (reference 8) that in the plane 
of the wing 


D 3 RD. ("ар а 


- (18) 


(2— z) ÊP - 
o o A ___ a 9) 
(y—y)N(z —2)—8((y—y) 


where the order of integration must be carried out as indi- 
cated (i. e., the integration with respect to yı must be made 
first). The letters R. P. mean that the real part of the term 
is to be taken. Such a device can be used since the double 
integral must always be a pure real quantity in the area т 
(Ар/9 is real everywhere on the plan form) and a pure 
imaginary quantity over the rest of the area indicated by 
the limits on the integrals (see fig. 4). The average vertical 
induced velocity along the span may be obtained by applying 
to equation (19) the operator 





1 
2тх |; 
and since 
Ар Awym E 
i= 20 
| @ үт? — y faa) м9 
equation (19) may be written in the form 
— c Вто fi), Wo f til, fi(zi)d zi 
Met 2 i dll 0 t— T] (21) 


fi 
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FiGURE 4.— Regions for which the double integral term in equation (19) із a pure гея] ог pure 
aginary. 


The derivation of J, is similar to that used for equation (9) 
and leads to the expression | 


neto _ а= Big ds 


I;— m R.P. BT oe IUE EMEN YOU ee EN 0 
Jaana) qv [mni 2) — т] (mG:— 2) 3-91] 


(22) 


It is possible to find an exact solution for f;(z) by means 
of equation (21), but the discussion of this analysis is reserved 
for & subsequent section. 

Rectangular pian form.—Equation (19) can also be used in 
the case of a rectangular plan form by an appropriate change 
in limits; thus, 


v Ар 
8 Fa— V. £ г 
uw — 14 PRP. dap dy; 
T 0 —t& 
Goce : 
g (23) 


(уу) 4(z — r} — Ву — yn) 


where again it should be stressed that the order of 
integration cannot be reversed. The regions in which 
the integration yields real or imaginary results are 
shown in figure 5. As in the case of the subsonic rectangular 
wing, the value of w will be obtained only along the 2 axis 
so that у in equation (23) can be set equal to zero. The 
loading will be assumed to have & form 


Ap 4 Ve гү EN [3 
"2 (а) y 





(24) 


which is similar to that used in the subsonic case except that 
the reference length is now the semispan instead of the chord. 
Such & difference is reasonable since in the supersonic case 
the position of the trailing edge cannot affect the loading 
'on the wing. 

Finally, therefore, when y—0 equation (23) becomes 


1058 . 





—Pure real 


Pure — —. 
imaginary 


/ 


ЕтасвЕ &.— Regions for which the double Integra] term fn equation (23) fs pure real or pure 


imaginary. 
Mero 
= ®\, Wo |? NSB/ 
=), (5) +55 | ——u dn (25) 
where J, is given by the equation 
1 m = 3 І 
L=R.P. f у dy, (26) 


The solution of equation (25) is deferred to a subsequent 
section. P 


LOADING ON WINGS 


The previous section was devoted to the development of 
the integral equations which are to be studied for the two 
types of plan forms in subsonic and supersonic flight. In 
order that this study can proceed in & natural manner, the 
arrangement of the presentation has been changed so that 
the plan form is the principal division and the speed is 
subsidiary. | 

TRIANGULAR WINGS 

Supersonic case.—The decision to solve for the loading 

on the supersonic, triangular, flat plate by analyzing equa- 


tion (21) was not an obvious one since the exact solution of _. 


the linearized partial differential equation for this case has 
already been obtained. (See, e. g., references 10, 11, and 
12.) Thus it is known before starting that the value of 
f;(z) in equation (21) must be 1/E where E is the complete 


elliptic integral of the second kind with modulus J1—m3g*. ^ 


However, these solutions were obtained by an entirely | 
different procedure so that by solving equation (21) and 
comparing the two results & check on the aecuracy of the 
method із obtained. Furthermore, when the subsonic 
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problem is analyzed the same general procedure will be 
followed and the results can then be accepted with greater 
confidence. 

The first step in the solution of equation (21), in which w 
has been set equal to w since the wing is a flat plate, is to 
change variables by the transformation &=2,/z. This gives 


8m ,, 1 (b foEdE 
t far | 2848 т, me 


In the equation for J}, the transformation my1=n/2 was used 
so that 


= n= (75. J(1— &y— Вт? т? 
1-6 т МЕЕ itn | 
whieh 18 completely independent of z. Тһе partial deriva- 
tive of both sides of equation (27) with respect to z gives 


AG) — 55 | 80088048 мы) 


Equation (28) is а homogeneous linear integral equation. 
The solution to equation (28) is simply f,'(z) —0 or, what is 
equivalent, fi(z) equals a constant, (fs). say. By means of 


(28) 


equation (27), this constant can be evaluated. Hence, 
TEM q 1 ГАТ алаа | 
(fs) =| =e +), = (29) 


which represents the solution to the problem. 
Г, was calculated analytically as in appendix A, and then 
the value of (}%)о, as given by equation (29), was determined 
by numerical integration. For Sm=0.8 the result of this 
computation was 0.708; whereas the true value given by 
1/Е is 0.705. 

Equation (27) can also be solved when the wing is slender 
with respect to the Mach cone by considering 8m to be 
small. Setting $8m—0 yields 


(+) dq 


"o "acd rU 


and this is readily evaluated to give 


(Laam о = a 


4y1—£ 


The integral equation reduces to 


.. (80) 


I= Nr hfdzt)éh (z£)dt 
E Jy1— £j 


which by a retransformation of variables z,—z£, becomes 


z= | лың | 
0 Jzi— xj 


The integral 


(81). 
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"Equation (31) is a special form of Abcl’s integral equation, 


the unique inversion * of which is, in this ease, f(z)—1. 
This is easily verified by direct substitution. 

The simplicity of this result is not accidental, of course, 
since the value of f(z) was originally introduced by equa- 
tion (20) as a correction factor to the slender-wing-theory 
solution. 

Subsonic case.—The study of the triangular wing pre- 
sented in the preceding section was made first at arbitrary 
supersonic Mach numbers and then at a Mach number equal 
to 1. In keeping with this order of decreasing speed, the 
subsonie flat plate will be studied first at sonie speed and 
then for general subsonic Mach numbers. 

An inspection of equations (9) and (22) is sufficient to 
show that (Л)вњ-о is equal to (Jz)emuo. Hence, equation 
(30) can be substituted into equation (10) and there results 
(since again 10 is set equal to wo) 


pue пл (dn nh (2 a) аа 


Ls ——_— ore 


rz Weca ж drar 


and this reduces immediately to 
21} 1 eA (2) а dx 
a 


fein? 

It is now obvious that equation (31), which was derived 
from supersonic wing theory, and equation (32), which was 
derived from subsonic wing theory, are identical. Clearly 
this establishes the continuity of the theory in passing from 
the supersonic to the subsonic regimes. 

The study of the general subsonic case leads eventually 
to the numerical solution of an integral equation. However, 
an idea of the qualitalive form which this solution must 
assume can be gained by some preliminary analysis. 

First write equation (10) in the form 


(82) 


(33) 





iuis f т2:} beata f ink (3) лема 


Ja — 21 То 13 


where zy —2jes A and ®/w=1. The evaluation of J, 
is given in appendix A, and a plot of 47, against za/zo for 
(8m)* equal to 0, 0.05, 0.10, and 0.20 is shown in figure 6. 
Obviously equation (33) is a singular integral equation. 
We have already seen that for 8m-—0 it is an Abel type 
integral equation with a ‘power singularity at the upper 


t If Abel's equation is written In the form 


fig)» - er = 


d 3M = 
* dj. Jz-y 


ita inversion ія 


£x) 


‘and this inversion ія unique for functions £(z) that fulfill the condition 


tim 
6-0 &(7-6) =й. а<>к<у (&) 


Since the solutions for the velocity throughout the flow field must also satisfy condition (2), 
the value of g(z) given above is unique in the class of functions available. 
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limit. For finite values of Вт, however, Г, is seen from 
figure 6 to be bounded and nonzero in the interval 
D <rt о; hence, for such Bm, equation (33) is a combi- 
nation of two types of singular integral equations, the Abel 
type and the Cauchy type, the latter having a first-power 
singularity in the interval of integration. 

Experience with the Cauchy type integral equation 


f= =| & n (34) 


which arises, for example, in the — of subsonic lifting- 
line theory and two-dimensional airfoil theory is useful in 
the present problem. Thus, the solution ta equation (34) 
is not unique (even when g(z) is restricted according to 
equation (а) in footnote 4) unless some additional condition 
is given. Such a condition might be the requirement that 
#(1)=0; this would correspond in two-dimensional airfoil 
theory to the specification of the Kutta condition at the 
trailing edge. Further, it is known that if the condition 
2017==0 is satisfied in equation (34), then g(x) tends to 
infinity as z approaches zero. 

Since the solution to equation (33) must also satisfy the 
Kutta condition, the above discussion leads one to anticipate 
for the shape of }, (20) а curve something like that shown in 
figure 7 (а). On the basis of such qualitative knowledge, a 


simple numerical procedure was set up and used to calculate 
This procedure, which is 


the solution to equation (33). 
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based on the assumption that у, (х0) is constant over each of 
nine equally spaced intervals, is presented in appendix B. _ 
The results of the analysis for (Вт)? equal to 0.1 (i. e., 
8.i=1.26) are shown in figure 7 (b). j 
In order to check the results derived from the method just 
mentioned, equation (33) has been solved in an alternative 
manner. In this second approach it is assumed that fi(z) 
can be approximated by a second degree polynomial which is 


multiplied by two factors, one that vanishes for z;,—1 and the _ 


other that tends to infinity as лу vahishes. 

First consider the behavior of А (х) for small x. 
study of the solutions to the two integral equations repre- 
sented by equations (32) and (34) will illustrate how ў (х0) 


А brief 


can be analyzed in the vicinity of the origin. First consider  . 


the Cauchy type integral equation (34) and assume that for x 
small g(r) сап be expressed in the form 

I as Я 
ЕП 24 nT, 


&(®)= 17>8`>0 


Then equation (34) becomes . { uo 





dz = L e*-idy 
2 a. f r—y 


z(y—z) #1 


E 
Ау) = |, 
which in the limit as y—40 reduces to 


р а, : 
Pel am Go 


t dr 
o z'y—z) 




















a 









| 


Figrre 6.— Variation of 1/2 with 1s/1y. 
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o 3 intervals 
А 6 intervals 
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e 
4 
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(a) Estimated, 
(b) Calculated. 


FIGURE 7.—Variation of f with ze. 


or, setting z—yt£, to 
lim 1" 
p y 0 


i (1—2) - fo 0*2: scs] 


The term on the left is indeterminate and implies 


lim (9 48 (7 
v9], PA- Je WI £) 


This equality is satisfied, since 5 is to be greater Шап zero and 
less than 1, only by the value 6—3$. The fact that this must 
be the exponent of 1/2 in the expression for g(x) сап be verified 
by inspecting the known solution to equation (34). 

A similar analysis applied to Abel’s integral equation in the 
form given by equation (32) shows for that case à must be 
zero which, again, agrees with the known inversion. — 

Finally, in appendix C this same approach is used to dis- 
cover the initial behavior of fı (zo) in equation (33). Figure 8 
presents the results of this analysis throughout the range of 
B.A for which 7, was calculated. 

Ав was already mentioned, the variation of file) in the 
vicinity of z,—1, that is, near the trailing edge of the wing, 
is fixed by the Kutta condition. А useful statement of this 
condition that holds in both subsonic &nd supersonic flow is 


un д==т cot [(1—)я] =0 


FIGURE 8.—Varlation of 5 with reduced aspect ratio А.А. 


that the magnitude of the loading at the trailing edge must 
not be infinite. However, the only pertinent solution to 
equation (34), and hence to equation (33), that is not infinite 
at z—1 is that which is identically zero there. Further, as 
the trailing edge of the triangular wing is more and more 
closely approached, it is reasonable to expect that the shape 
of the load distribution in its vicinity approaches that in the 
vicinity of a simple two-dimensional wing trailing edge, the of- 
fect of the wing plan form vanishing as the ratio of the dis- 
tance from the leading to trailing edge tends to infinity. On 
the basis of these considerations, fi(zo9) should approach zero 
as the term (1—2,)% approaches zero when хо tends to 1. 

It is also apparent, however, that f(z) equals 1 for all 2, 
between zero and 1 when the Mach number is unity. Fur- 
ther, as Mae approaches 1 or as the angle of sweep approaches 
90°, the effect of the plan form on the shape of the loading 
near the trailing edge becomes increasingly important or, in 
other words, the trailing edge of the triangular wing must be 
more and more closely approached before the shape of the 
two-dimensional load distribution is simulated. An ex- 
ponent to (1—zx,) which satisfies this requirement as well as 
those in the preceding paragraph is y/2(y 4-1—24) where y is a 
function of Mach number and vanishes as M41. 

Finally, therefore, it is assumed that f (xo) can be expressed 
by the equation 


filz) = 


(ауана 
ry 


(Go ttot a28) (35) 


The values of the constants do, a1, €a, and y are determined by 
satisfying equation (33) at four chordwise stations and are 
given for (Am)? equal to 0, 0.1, and 0.2 in table I. The 
accuracy to which these solutions satisfy the integral equation 
is indicated by table II. 

A comparison, for (8m)*=0.1, between the solution given 
by equation (35) and that derived by the method outlined in 
appendix B is shown in figure 7 (b). "The strictly numerical 
method presented in the appendix was used for three different 
interval spacings, results for which are indicated in the fig- 
ure by the symbols. Presumably the accuracy of the method 
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increases with increasing number of intervals used. А study 
of the figure shows that the numerical method is apparently 
approaching * the solution given by equation (35). Subse- 
quent values in this report are based on the solutions repre- 
sented by the latter equation. 

The final curves for fi(z;) are shown in figure 9. А discus- 
sion of the integrated values of the loading will be given later. 


RECTANGULAE WINGS 


The discussion of the triangular wing was divided accord- 
ing to the Mach number. "The same division will be used in 
this section, starting with the discussion of the results for 
supersonic speeds, then with that for both supersonic and 
subsonic theories at sonic speeds, and closing with a discus- 
sion of the subsonic development. 





FURE 9.—Final values of fi(ze) given by equation (35) and table I. 


Supersonic case.-—The solution of equation (25) will give 
the loading on а rectangular wing flying at а supersonic 
Mach number. The evaluation of the integral J, is carried 
out in appendix A where it is shown that J, can be expressed 
in terms of complete elliptic integrals of the first and second 
kinds. Having the expression for 7,, а numerical solution 
may be obtained for f,(z/«8) (see appendix B). Figure 10 
shows & plot of f, (& factor representative of the chord lift 
distribution) for a flat plate wing as a function of 2/sf, 
the ratio of the distance back from the leading edge to the 
magnitude of the reduced semispan. The value of f, given 
by equation (25) can be checked in the interval 0 € (z/s8) <2 
because the exact solution to the complete linearized partial 
differential equation can be readily obtained there. The 
comparison is given in figure 10. The fairly rough agreement 
shown is not surprising since equation (25) is derived on the 
assumption that the spanwise distribution of load is elliptical 
at every chord station, and certainly this assumption is 
least accurate in the interval where the comparison with 
the exact results is made. The area under the exact and 
approximate curves in figure 10, between the initial value 
and that at which f,—0, is nearly the same. (See the next 


+ The maximum difference in. total lift is already less than 6 percent (based on the lower 
value) and the difference in center of pressure Їз negligible. 


272483—54——688 
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section on Aerodynamic Characteristics.) The integrated 

value of f, as given by equation (25), therefore, can be used 

for z/sB>2. | 
As for the qualitative nature of the variation, figure 10 


shows that the loading on & narrow rectangular wing flying 
at supersonic speeds falls linearly to zero, becomes negative, — 


and then oscillates between negative and positive values, 
the amplitude of the oscillation being so heavily damped 
that after the third change in sign the magnitude 1з practically 
zero. 

It should be noticed in studying the results of figure 10 
that the entjre resultant lift of the wing is concentrated 


— ——— MH 


in the interval 0< (2/88) X2. But as the Mach number. .... 


approaches 1 this interval approaches zero, and the entire 
lift of the wing is carried in a strip along the leading edge.’ 
Such а solution violates, in the vicinity of the leading edge, 
the assumption on which the theory is based and should 
be considered only as a theoretical limit. 


Results for the lift and pitching moment on the rectangular _ 


wing will be developed in a later section. 

Subsonic case.—The study of the subsonic rectangular 
wing stems from equation (16). The first step in the analysis 
of the equation will be to consider its solution at fs=0 
and show that this is continous with the supersonic results 
there. 

The value of J; can be written (equation (15)) as 


(^ Ка х1) 82у? ш 
ne fy tage % 


and for Ss=0 this becomes 


(z—aiye 
һ=| 


(х—т)т 


Tı «zr 


21 2> 2 


and hence equation (16) сап be written 


Wa [F 2 
w= һ(& dii 
Co Jo C 









UE 
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Exact linearized 
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FiauRE I0.— Variation of chord wise correction factor fi for supersonie rectangular Wing. 
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t This result also follows by inspecting equat fon (25) for the values £10. 
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Equation (86) is identical with the form of the supersonic 
equation (25) at 8s—0 so that once again the continuity 
of the subsonic and supersonic theories at the sonic speed 
range is established. Furthermore, equation (36) shows that 
if w/w, is constant then (20) must be. zero everywhere 
except at points where it can be represented by ә pulse, 
the integral of which has a finite magnitude. From the 
supersonic discussion, it is clear that one such pulse exists 
and is located at the leading edge. 

The evaluation of 7; for Bs>0 is given in appendix A. 
The numerical solution to equation (16), assuming the Kutta 
condition at the trailing edge, is given in appendix B for 
values of reduced aspect ratio BA equal to 0.33, 1.0, 1.5, 
and 2.0 For an aspect ratio equal to 2, these values corre- 
spond to Mach numbers of 0.986, 0.866, 0.662, and 0, respec- 
tively. The results of the computations are shown in figure 
11 where the chord wise lift distribution factor ў (20) is plotted 
against zo for the various values of pA. Ву comparison of 
figure 11 with figure 10, it can be seen that in the subsonic 
case the loading drops monotonically from infinity at the 
leading edge to zero at the trailing edge and does not oscil- 
late in the afterportion, as in the case of the supersonic wing. 

When 8 equals one, these results can be compared with 
those obtained by Wieghardt and presented in reference 9. 
Figure 12 shows the comparison for two values of the aspect 
ratio. Curves are also shown in the figure for the loading 
obtained by using the method given in appendix B but by 
satisfying the integral equation at only six and three points. 
The latter curve is in better agreement with Wieghardt’s 
result and, since Wieghardt (although using a different. 
method involving Birnbaum functions) used only four points, 
this may account for the discrepancy between the final 
results of this report and those of Wieghardt. 


2.0 
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FiauRx H.—Variation of chordwise correction factor fs for subsonle rectangular wing. 
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Wieghardt 


This report 
final 





Fringe 12.— Load distribution along center line of rectangular wing. 


AERODYNAMIC CHARACTERISTICS 


The previous section presented solutions for the loading 
on triangular and rectangular wings flying at subsonic and 
supersonic speeds. This section will be devoted to the con- 
version of these loadings to expressions for lift and center 
of pressure. 

LIFT 


By definition the lift coefficient can be written 


1 А 
c.g [|22 dxdy 


and this will be evaluated for the various cases for which 
the loading coefficient hes been obtained. 

Supersonic triangular wing.-—Since the exact linearized 
value for the loading on the triangular wing flying at super- 
sonic speeds has been derived, the lift coefficient. can be 
written in the form 


(37) 


С ; 
on es 


where A is the aspect ratio and E is the elliptic integral of © 
the second kind with modulus k= 4/1— 8*m*. 
. Subsonic triangular wing.—In the case of the subsonic 


triangular wing, equation (37) becomes 
T 
" 1 (2) ау 


(^а [i fomir 


and this ое (since A=4m) 





й 
б, [ Zo f i(zo)d zo (3 9) 


The numerical evaluation of equation (39) is not difficult 
since zo fi(zo) vanishes at z9—0. 

Supersonic rectangular wing. --For values of eq« 248 the 
exact value of the lift coefficient on a rectangular wing flying 
at supersonic speeds has been obtained and ean be written 
in the form for 84751 | 


C, 4pA— | 
«А. BAY 40) 
When BA<1 equation (37} must be used in connection with 
equation (24) and there results 
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LE E © Я i Gg Pov а 


which reduces to 
CEN fs E ы. 


and this сап be written In tbe form 


xr E fe (вва) ee 


57 


m LL 
which becomes, if 2 = 82 


for BA<1 


=> | nest (41) 


Subsonic rectangular wing.—The equation for the loading 
on a subsonic rectangular wing, equation (14), placed in the 
formula for lift coefficient yields 


г ' fa Ж Er nue 
=з], 4 2 ia (22v л 


which becomes 


6. 2 JN fa(xo)d 1o (42) 


The evaluation of equation (42) by numerical means re- 
quires special consideration since (х) approaches infinity 
at the leading edge as shown in figure 11. To this end, re- 
write equation (42) in the form 

C zf" m [1 

i5], f Ка) |, Ро) zo (43) 
and equation (B4) in the appendix (for the special case in 
which z421) in the form 


== f ма) +1805 Es drgd- 


+ an fale) | r+ | dz, (44) 


An application of the mean-value theorem yields 


|: fx(xo)d zo 


Tag {зк— f Га) RE gom E dzo} (45) 
8 } 


where E, has the modulus $s which а В414/4(1-#)-+ (В.А)? 
und where 0<0<є. The combination of equations (43) 
and (45) y1elds an expression for the lift coefficient involving 
only the load distribution from а distance e/cg back of the 
leading edge to the trailing edge. 





PITCHING MOMENT 


By definition the pitehing-moment coefficient &bout the 
upex or leading edge and based on the root chord can be 
written 


CORRECTIONS TO SLENDER-WING THEORY 


(46) 


Equation (46) will be applied to the various loadings which 
have been studied. 

Supersonic triangular wing.— he exact heared value 
for the pitching-moment coefficient on a triangular wing 
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flying at supersonic speeds has been derived elsewhere and - 


can be written in the form 


Cn 


psi (47) 


NM n 
3E 
Subsonic triangular wing.—The derivation of the pitching 
moment on а subsonic triangular wing proceeds in the same 
manner as the derivation of lift and there results 


Cn = rf Lo fi(zg)d za 


(48) 

This expression can be easily integrated numerically. 
Supersonic rectangular wing.—For values of 8.4 greater 

than 1 the pitching-moment coefficient on a rectangular wing 


is given by the equation, for BA>1 


Cs. 68.1—4 
zd Am * (49) 


When &A«1 the solution to the integral equation must be 
used and the final expression can be written 


he [ч zof(x2)d xa 


Subsonic rectangular wing.— The equation for the pitching- 
moment coefficient on а subsonic rectangular wing follows 
in the same manner as did that for the lift coefficient. 
Henee, 


Cn 


ad (50) 


( 1 
a-l, 105020) zo (51) 
and, since the variation of хо fa(zo) is as indicated in figure 13 
the numerical integration of equation (51) is simple. 


CENTER OF PRESSURE 


Since the pitching moment is based on the root chord, the 
center of pressure of all wing plan forms can be written 


EM, Ca 
Со 1, 


(52) 


DISCUSSION OF RESULTS 


Figures 14 and 15 show the variation of the lift coefficient 
and center of pressure on triangular and rectangualr wings 
for values of 84 between zero and 2. For the triangular wing, 
the differences between the subsonic апа supersonic cases 
are not large in this interval of reduced aspect ratio; the 
subsonic wing develops somewhat less lift and its center of 
pressure moves forward as В.А increases. The characteristics 
of the rectangular wing, however, show a large variation in 
passing through the speed of sound. 

The subsonic rectangular wing has a variation of Cz/da 


with 48 which is almost identical with that for the subsonic 
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Reduced aspect ratto, A 


(a) Lift. 
(b) Center of pressure, 
Fiaure 14,—Aerodynamic characteristics of a triangular wing having в low value of 8.4. 


triangular wing. Unlike the triangular wing, however, tho 
curve for z,, р./со on the rectangular wing shows this lift to be 


carried farther and farther forward with decreasing 84, from. 


the quarter-chord position at #А= œ all the way to the 
leading edge at BA=0. 

As the speed is further increased and the АЕ wing 
enters the supersonic speed range, the magnitude of the lift 
begins to oscillate with increasing amplitude. This continues 
until the reduced aspect ratio rises to one, after which, as 
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(b) Center of pressure. 
FInuRE 15.—Acrodynamic characteristics of rectangular wing having a low value of 8.1. 


В.А increases still farther, the value of Cz/Àa falls uniformly 
to zero according to the expression 


CG, 4 1 
=F (17333 2 i) 8A21 (53) 


which is the exact equation for the lift coefficient given by 
linearized lifting-surface theory. The difference between 


- the values of lift coefficient given by equations (41) and 


(53), represented in figure 15 by the solid and dashed lines, 
respectively, has already been discussed in the section on 
loading; the approximate solution is based on the assumption 


. that the span loading at each chord station is elliptical and 


such an assumption is increasingly unrealistic for increasing 
BA. The dotted line shown in the figure appears to be a 
reasonable interpolation between BA=1, the lower limit to 
which equation (53) applies, and a nint where the approxi- 
mate solution given by equation (41) сал be considered 
trustworthy. 

The variation of the center of pressure on & supersonie 
rectangular wing indicates that the wing is unstable for all 
positions of the pivot point behind the leading edge for 
values of B84 around 0.4, the center of pressure, in such а 
range, having moved forward of the wing leading edge. 
Аз 84 increases past the value of 0.5, however, the center 
of pressure moves back along the wing and rapidly ap- 
proaches the midchord point, its location for a wing of 
infinite aspect ratio. 


L 
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COMPARISON OF RESULTS FOR SUBSONIC TRIANGULAR 
WINGS WITH OTHER THEORIES AND SOME EXPERIMENTS 


DISCUSSION OF THE THEORETICAL RESULTS 


Several other published theories can be used to calculate 
the forces апа moments on low-aspect-ratio triangular wings 
flying at subsonic speeds. А comparison between values of 
Ca, and center of pressure given by those theories and by 
the method of this report is summarized in figure 16. 
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Етогих 16.— Various theoretical and experimental results for triangular wings showing the 
variation of BCL, and ex feo with reduced aspect ratio, A.A. 
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The theories of Falkner (reference 13} and Weissinger 
(reference 14) are well known and will not be discussed in 
detail here. 
were presented by De. Young and Harper in reference 15 
for aspect ratios equal to 1.5, 2.5, and 3.5; those obtained 


The results obtained from Weissinger’s method 


LUG] 


from Falkner’s method were presented by Berndt in reference . 


16 for aspect ratios equal to 1.0, 1.67, &nd 2.5, and by 
Falkner in reference 13 for an aspect ratio equal to 4. Both 


Berndt and Falkner used 126 vortices and 6 control points. 


so the accuracy of their calculations is similar. It should 
be mentioned, however, that in order to make Berndt's 
results consistent with Falkner’s, the values of С. given 
by Berndt have been multiplied by a factor " suggested and 
used by Falkner in reference 13. 

The results shown for Lawrence's theory were presented 


in reference 17 and are based on а method originaliy given | 


in reference 18. 


The differences between values obtained from each of the 
theories can be attributed to differences in the various simpli- 


fying assumptions used. The Weissinger and Lawrence - 


theories agree on values of Cz_ but disagree as to the position 
e pc cc NN 


of the center of pressure. Faikner's method and the method 


of this report yield results in good agreement for both Cz. - | 


and center-of-pressure location. However, for 8.1—2, the 
Weissinger-Lawrence value of Cz, is about 9 percent below 
and the center-of-pressure position about 3 percent (based 
on the root chord) farther forward than similar values. 
obtained in this report ê or by Falkner's method. 


DISCUSSION OF THE EXPERIMENTAL DATA 


Experimental data for low-aspect-ratio triangular wings Ы 


are given in references 13 and 19 through 22. "The sections, 
aspect-ratio range, Reynolds numbers (based on the mean - 
aerodynamic chord), and section thickness ratios are piven 
in table III. 
The source of the experimental values for the d=4 wing 
(given in reference 13) is some unpublished British wind- 
tunnel data. The Reynolds number is given as “high” and 


the section is not specified. Reference 21 presents the results — 
of experiments made in the Langley free-flight tunnel on . 


some flat-plate (% inch thick with rounded leading edge) 
models having beveled trailmg edges. Hence, the values for 


_the thickness ratios listed for these tests are effective thick- 


ness ratios, being, in fact, the plate thickness divided by the / 


root mean chord of the wing, (2/8)¢o. 


The experimental date presented in table IV are portions — 


of the data given in references 19 and 20 and the values 


shown in figure 16 were obtained from curves constructed ~~ 


by means of these numbers. Reference 21 presents experi- 


mental results in graphical form and the slopes used in 
figure 16 were read from these graphs. It should be men- _ 


tioned that reference 21 also gives results for ап aspect 
r&tio 0.5 wing but the data presented are not sufficient to 


т According to Falkner (reference 13): “It has not been possible to establish the factor for ` 


all cases, but figures derived from a deita and other wings suggest that the factor can be taken 


tentatively as independent of aspect ratio, and to vary 8s 170.029 (tangent of sweepback of © 


quarter chord)." 

E Tt can be shown, however, that the Integral equation used by Lawrence to obtain the 
curves shown in figure 16 and the one used in this report both approach the same values of 
er. and z«,, [с as the angle of sweep goes to efther x/2 ог zero. In the formar case they approach 


the familiar Jones" low-aspect-ratio-theory results, and In the latter case they yield 3r for Cr. 


and 1/3 for re. Jeu. 
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fix the slopes of the lift and moment curves near zero angle 
of attack. Hence, no values for this aspect ratio are pre- 
sented. Finally, the values taken from reference 22 were 
read from graphs presented therein and based on integrated 
pressure distributions (measurements taken along five span 
stations). 


COMPARISON OF EXPERIMENT WITH THEORY -LIFT-CURVE SLOPE 


The experimental data shown in figure 16 (а) do not seem 
at first glance to favor either group of theories. However, 
since all these theories are based on the assumption of zero 
thickness, 16 1s pertinent to examine the data on the basis 
of thickness ratio. A correlation on this basis is presented 
in figure 17. It shows that with decreasing thickness ratio, 
the ratio of the experimental values to those predicted by 
the method of this report (in the range 08431 x2) or by 
Falkner's method (in the range 2844) tends to unity. 
Of course, this trend is not conclusively borno out by these 
comparisons and should be tested by more experimental 
measurements in the lower thickness ratio range. 

It is interesting to notice that a good approximation to 
the results caleulated in this report and by Falkner’s method 
is given by the equation 


Cp 2rcA _ | 





(54) 


pAt+2 





о .02 0 4 06 08 ло 12 14 
Thickness ratio, 7c 


Fiavre 17,— Varlation of (Cz) „(К Cz, ) ,, , the ratio of the measured lift-ourve slope to that 
given by the theorot{ca] methods of Falkner or this report with thickness ratio. 
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where 


NULL semiperimeter _ | EIU 


- wing span 


Equation (54) was derived by R. T. Jones (see referenco 23) 
as a first-order correction to the value of Cz, given by 
lifting-line theory for wings having elliptic plan forms. It 
has been found, however, to provide а good estimate for the 
aspect-ratio correction to wings of various plan forms. For 
the particular case of the triangular wing in compressible 
flow, equation (54) becomes 


BC, = кВА 


1484.) 1 (82у 


and its variation for 0<8Д <4 is shown in figure 16 (a). 





(55) 


COMTATISSUN OF EXPERIMENT WITH THEORY-CENTER-OF.PRESSURE 
LOCATION 

The comparison between experimental and theoretical 
values of center-of-pressure location is shown in figure 16 (b). 
In weighting the experimental points shown, it should be 
remembered that the tests recorded in reference 21 were 
made on a flat plate with a beveled trailing edge. When 
such a wing is at a positive angle of attack, the bevel causes 
the flow to separate at the surface discontinuity producing 
an effective upwardly deflected flap. The result is a loss of 
lift on the rearward portion of the wing and a forward shift 
in center of pressure. Hence, the experimental values for 
Z5.» ./€o taken from reference 21 should be low relative to the 
values given in the other references. With this taken into 
consideration, the experimental values shown in figure 16 (b) 
are in fairly good agreement and again favor, at least, the 
theoretical results based on Falkner’s method. 


Ames AERONAUTICAL LABORATORY 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
Morrerr Fieup, Cauir., Nov. 28, 1952. 


* This approximation із suggested in an article by R. T. Jones and Doris Cohen, "Acrody- 
namics of Wings at High Speed," to be published by the Princeton University Press in 
section A of the book entitled Applied High-Speed Aerodynamics, volume YI of the High- 
Speed Aerodynamics and Jet Propulsion Series. 


APPENDIX А | 
EVALUATION OF SPECIAL INTEGRALS 


THE INTEGRAL 7, 


The evaluation of J, will be discussed first for the case in 
which œz, and second for the case in which z«z,. 


CASE I, z»a 
It is possible to write J, in the form 


eere r a = =——зиш=————шш=——=ч 


E = Al 
Jo qy (n шо) (Qu — 9) x 


where gy m(x—2,) and gu —m(-FEz). The linear term in 
the lower radical of the integrand can be eliminated by the 


transformation n= (еє-Е8ї)/(1-Еї), and the integral becomes 


в 





ó—a 


eS s 
(рое) (e — ду) J 5-50 
Боё 


(ио a383 m3) + (p2+ ep m?) 


а ЕЕ ар (А2) 
росе = на 
(14-0 3-041 3 ( = G=) Р 
where 
LL Bano bmm) + uo v U + Вт?) (ud Вт 2) 

m | B*m*(u d ug | (A3) 
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з но B^m!ui)— tox 1 3- Bm?) (ud Вт n) (Ad) 


B*m*(ui4- н) 


The expression for с and ё may be combined to give the 
useful identities 
ug = — esp m? 


(ш— €) (8— pr} + (1 — 0) (— ро) = 0 


Using fundamental properties of even and odd functions, 
equation (A2) may be reduced to the form 


- Bo T c?f^m! E 1 fo 1 
Bac cmm 


1 

tap 
ЖИ; i 
y Г | 


by the substitution 
рі ô 
bio 
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u= 








and where 
B (ши є) | 
(ui? — eò) (8— v) 


By introducing the Jacobian elliptic functions in the 
transformations w=enu, the integral reduces to | 


Е. < 
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(А5) 
where 
т ` 
20)" 
PL = = ( (Аб) 
к= i 
and 
ron 
pi 0 (АТ) 
E k? 
The integration may now be completed and 
T _2 | po б6°8°т E [уу +Ё Gp a)- 
: kN (u1— a) (с — go) A, l+» NINE 
с vith Р» 
PV TE O (от) | 
_ Y УІ Ho LZ 
ЕЕЗ 
where | 
G(x)-- EF" (x) +R’ (x) КЕ (x) (49) 


the modulus of the elliptic integrals being ог k. 


CASE 2, x€x 


The procedure for obtaining the solution for J; in this case 
їз identical to that followed in case 1 except that in order to 
fulfill the contition that >—1, с and & must be defined in 
the following manner | 


вора Bm). пау ЕВ) (a E Ветар), 
$ = Вт? i ро) - (410) 
раро Bm?) — Hoy (1 + 83m?) (pd +t Вт? шї) pm 
i, (00 Вто) Аш _ 
In this case it can be shown that in equation (А5) »,750 


and м< —1, and the solution for Г, is 
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where G(z) is defined as in equation (А9). 
THE INTEGRAL I, 
Writing J, in the form 
G aF WR. (А18) 
52—31° 
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сти, k= ES E 
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Г, can be integrated to give 
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reduces equation (А15) to the form 
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The integration may now be completed and 
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In this case 
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and the transformation 


is made where 








Equation (A19) then becomes 


ае: 39 
боз 
The integration can be M, so that 


kcal &K(:— EB -1) 


sacar eon 
fe 28(1—68)») 1—8,4-28,:3/ 


1—6,4-28,k* С 3465  ——— 


bo | 201—0) iS =), 
DV V xz } 





— 
a (588) A20 
THE INTEGRAL I, 
CASE 1, 0€ Xi X —fs 
The integral Z, can be written 
53,3 | B 
[ез E goa re a2 umi Jo Ya =F йу (А21) 


When the transformation Pm is made, the expre ssion 
may be integrated to give 


Теа) |. dn?udu-2(x—az))f (A22) 


where 
828? 
2, 
(2—21) 
CASE 2, X-fs« Xi€ X 
In this case J, сап be written 


1—1 З 
v <z By? 
т=з Í с 208 dyi 


The transformation enu-—fy/(r—2z)) applied to equation 
(A24) yields 


(A23) 


(A24) 





Ен ү 
=з | XO en'udu-s 2002—00) В, (A25) 
where 
bm (A20) 
and 
B. ES -(A27) 
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APPENDIX B 
NUMERICAL SOLUTION OF INTEGRAL EQUATIONS 


SUBSONIC TRIANGULAR WING 


Since the integral Л is a function only of the ratio ау, 
equation (10) can be written for w= 1 


1 == 1 Ышы уз E 


2 0 4i—g 


Е i Ee)dt (ВІ) 


where f=2,/x. It is now assumed that fi(x£/eg) may be con- 
sidered constant over small intervals. This reduces the 
solution of the integral equation to the elementary problem 
of solving a system of simultaneous algebraic equations. 
On the basis of such an assumption, equation (B1) becomes 


1 n 
l—» £D 


where 





ч 
2i—1N (77i © Е 
= Iaa 208948; j—1,2,8,..m (B2) 
udi ern d 














E he; 05:91 | 


Co 
«ri«— 
1<#< af 


The function 2,70) was calculated by numerical integra- 
tion and systems of simultaneous equations were obtained 
for values of n equal to 3,6, and 9. Solutions were found 
using the Gauss-Seidel method (for which the simultaneous 
equations were well suited). 


| (83) 
na h; 


SUBSONIC RECTANGULAR WING 


Substituting the value of J; given by equation (А14) into 
equation (16), one has for w=, 


(B4) 


A CERO 
where 
(EE art a 
and where l 
88 BA 


сату (22) вт» 


A satisfactory numerical solution of equation (B4) requires 


the solution of the system of simultaneous equations of the: 


form 


aba Quit). 

















AGA 


(B5) 


АС) m a) tC 


љав (25521) 


The convergence of the solutions to equation (B5) is indicated 
in figure 18 where"the value of n was successively taken to be 


3, 6, and 9. 


О 3 Intervals 
4 6 intervals 
g9 Intervals 





Етеган 18.—Solutions obtained for the variation of f; with zv using 3, 6, and 9 intervals tn B 


equation (B5). 
SUPERSONIC RECTANGULAR WING 


Equation (25) can be written when w-wy as 














l- f, SH £ (= =) (2 36 A (B6) 
where from equations (А22) and (А25) ' 
І — 2 (Pe F A | 
«(5 sa) se 4 (ВТ) 
2 pg n a = <i 
and where 
8B . „ХХ: 
imme. kz 8В 


By the application of the trapezoidal rule for. numerical 
integration to equation (B6), it is possible to write Fi(z/sB) 
explicitly as 


8) оон) 89] те 


where A(z/Be) is the interval of the trapezoid. 
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APPENDIX СЄ 


EVALUATION OF $ 
Consider the integral equation 


26m |" TT fi (2а) d25 | 


42 2—25 o Le 


In order to study fi(z;) in the neighborhood of the origin, 
first set r=% во that equation (C1) becomes 


ош f. r£ fi (xot) d£ Иш EI (B fi (tok) dE 
(1-8 A 1—£ 


for small values of zs, filza) can be 


(C2) 


then assume that, 
expressed as 


Go . EE 
ho) ту E E 2a аел. Е 


From this one has 


-IdE 4 ио EST (E) dE 
2r=lim | x f b ou €» 
which immediately implies 
| -idg | (^ Ends _ 
Tj oe pm == (C4) 


The value of à in the interval O¢S8A<2 that satisfies 
equation (C4) was determined by numerical integration. 
The results are shown in figure 8. 
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TABLE I.—YALUES OF CONSTANTS FOR EQUATION (35) TABLEIV.—EXPERIMENTAL DATA ON LOW-ASPECT-RATIO 

: TRIANGULAR WINGS TAKEN FROM REPORTS. BY 
BERNDT (REFERENCE 19 AND LANGE AND WACKE 
(REFERENCE 20) (SEE TABLE III) 


(в) Reference 19 (Ca referred to іса and taken about point ics back from apex) 


- 0—— M — mmi a — 


А = 503 | А. 5/3 Ami 















TABLE IL—VALUES CALCULATED FOR RIGHT SIDE OF = ш 0. 029 -4 8 A 165 0.036 - 0. 097 0.028 
EQUATION (33) USING EQUATION (35) AND TABLE I FOR ET |: 055 ou | -28 | —.095 1022 | —23 | — 09 . 012. 
f, (же) —.6 | —.07 —LT | —.0@ „015 | -L2 | —.08 „007 

[Esact solution would, In each case, yteld 1.000] S c i А6 E 1 E Um os Ex uc ine 
2.5 .lt5 | —.017 L3 043 | —.005 1.8 049 | —.009 
3.5 .1l60 | —.024 1.3 078 | —.012 2.8 064 | —.014 
4.6 .305 | —.081 14 114 | — 009 3.8 — 
| 4.4 149 | —.026 4.8 110 | —.025 


єз 


(b) Reference 20 (Ca referred to {се and taken about point fes back from apex) 


LI 
„2 
.4 
4 
‚5 
-8 
.8 
.9 
„95 





TABLE IIL—RANGES OF PERTINENT PARAMETERS IN 
EXPERIMENTS 










Reference 
number 
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(0 


— ка 
«aq 




















Swedish 3x10! 
: FRA (04-5106 
20 МАСА 0012 1.7 to 3X10 
21 (1) -4 to 9X 10€ 
22 МАСА 0012 2. 4X 10€ —À 


т See text for explanarion or qualification. 


